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Abstract. With the linear sigma model, we have studied Bose-Einstein condensation 
and the chiral phase transition in the chiral limit for an interacting pion system. A 
fi — T phase diagram including these two phenomena is presented. It is found that the 
phase plane has been divided into three areas: the Bose-Einstein condensation area, 
the chiral symmetry broken phase area and the chiral symmetry restored phase area. 
Bose-Einstein condensation can happen either from the chiral symmetry broken phase 
or from the restored phase. We show that the onset of the chiral phase transition is 
restricted in the area where there is no Bose-Einstein condensation. 
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1. Introduction 

In relativistic heavy ion collisions, large amount of pions can be created in small space- 
time regions. Theoretically, it is possible for these pions to condensate into the zero 
momentum and form Bose-Einstein condensation(BEC) at certain temperatures and 
densities. For this problem, there are some discussions in the literature, for example, 
the BEC of pions in the non-equilibrium process PjJ |2] , the pion condensation in thermal 
equilibrium incorporating boundary effects [HJ, and the search for the "cold spots" in 
multipion system |U|3], etc. These studies show a great interest on this issue. On the 
other hand, in high energy physics, the chiral symmetry which is spontaneously broken at 
low temperature is expected to be restored at high temperature and/or density j^j. Long 
before, Kapusta and Haber had discussed the relation between BEC and the spontaneous 
symmetry breaking by using different scalar models Their main interest, however, 

was focused on the physics of the electroweak system. Here in the strong interacting pion 
system, one may wonder whether there are some relations between the chiral symmetry 
restoration and BEC. 

Our motivation in this paper is to discuss the relation between these two phenomena 
in the pion system and find the phase diagram of them through a model study. As 
far as we know, this topic has not been much discussed in the literature. In 
the chiral symmetry restoration is discussed in the interaction pion gas through virial 
expansion, where BEC is only discussed in the case of free pion gas. Here we will 
discuss BEC and the chiral phase transition in the interacting pion system with the 
linear sigma model. The linear sigma model has been extensively used to discuss the 
chiral symmetry restoration at finite temperature fUJ El E] In our discussion, to 
address BEC and the chiral phase transition simultaneously, we will consider a finite 
pion density as well as the temperature. The pion chemical potential can be introduced 
in different ways. Theoretically, at given conditions when the pion number can be 
regarded as being fixed [12], one could introduce the pion number density. However, 
generally speaking, pion number is not a conserved quantity and one should find a 
certain conserved charge of the system. For the linear sigma model which has the 
isospin invariance, one can introduce a chemical potential associated with the isospin 
charge. For example, in [14] . the pion isospin chemical potential is introduced by the 
third component of the isospin charge and applied to the discussions of the thermal 
pion mass within the chiral perturbation theory. In our discussion, the isospin chemical 
potential of pions is introduced to the linear sigma model in the same way. 

The organization of the present paper is as follows. In section 2 the effective 
potential for the linear sigma model is calculated within the Cornwall, Jackiw, Tomboulis 
resummation scheme at finite temperature and finite pion density in the Hartree 
approximation. The gap equations are derived by the stationary condition, which form 
the basis of our numerical calculation. In section 3, the BEC and the chiral phase 
transition are studied by a numerical calculation. We discuss their [i — T phase diagrams 
and give the physical analysis. The last section comprises a summary of our results. 
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2. CJT formalism and effective Lagrangian 

In studies of phase transitions the effective potential is an important and popular 
theoretical tool. It is defined through an effective action which is the generating 
functional of the one particle irreducible graphs. A generalized version is the 
effective potential for composite operators introduced by Cornwall, Jackiw and 
Tomboulis(CJT) [12]. This formalism assumes that the effective potential V(<j),G) 
depends not only on <f>(x), a possible expectation value of the quantum field but 
also on G(x,y), a possible expectation value of the time-ordered product T<&(x)$(y). 
Physical solutions demand minimization of the effective potential with respect to both 
4> and G, which means 

dVj^G) dV{4>, G) _ 

d<p ' dG { ' 

This formalism was originally written at zero temperature. Then it has been extended to 
finite temperature occasion by Amelino-Camelia and Pi for investigations of the effective 
potential of the A0 4 theory [T^]. It has been also used in the linear sigma model for 
discussing the chiral phase transition at finite temperature At finite density as well 
as finite temperature, we will employ this formalism to discuss both BEC and the chiral 
phase transition with the linear sigma model. 

For the linear sigma model, here we will ignore the fermion sector. The Lagrangian 
can be written as 

(2) 

where a and tt are the sigma field and the three pion fields (71*1, ^2,^3) respectively. The 
explicit chiral symmetry breaking term is e = j^rn^^ where f n = 93MeV is the pion 
decay constant. The coupling constant A and negative mass parameter m? of the model 
are chosen to be A = 3(m^. — m 2 )// 2 . and — m 2 = (m 2 — 3m^)/2 > 0. In our discussion, 
we only consider the case in the chiral limit e = 0, then at zero temperature the pion 
mass is = and the sigma mass is taken as m a « 600MeV. 

The conventional treatment of spontaneous chiral symmetry breaking of the linear 
sigma model is to shift the sigma field as o — > a + 0, where do is expectation value of the 
sigma field and the order parameter of the chiral phase transition. Then the Lagrangian 
results as 

C = -(On) 2 + -(da) 2 - -(m 2 + -<j) 2 )a 2 - -(m 2 + -0 2 )tt 2 

_ I m 2 2 _ _ A( a 2 + ^2 )2 _ A 3 _ A 2 

2 y 2V 24 v J 6 r 6 r K ' 

In the following discussion, we will adopt Hartree approximation which means only the 
bubble diagrams with four-point vertex need to be calculated. So the last two terms in 
the Lagrangian are neglected |lUj . 
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The chemical potential is introduced through the third component of isospin charge. 
The Lagrangian (3) is invariant under the isospin transformation, therefore according 
to Noether's theorem the conserved current is = tt x d^fc and the conserved isospin 
charge is Vq = 7? x <9o7r. Then the partition function of the system can be written as 



[dUtf] [dU a ] [dn] [da]exp 



d 4 x(C + jlj • Vq 



(4) 



where LT^ and Tl^ are the conjugate momenta of 7? and a. The isospin chemical potential 
fli = (0,0,//), in which /i is the third component of isospin chemical potential, is 
introduced associated with V [141, thus 



jl! ■ V Q = /i(7ri<9 7r 2 - 7r 2 9 7r^ 



(5) 



Here it is more convenient to work with fields possessing well-defined charges: 

7T = 7T 3 . (6) 



7T_ 



1 , ■ s 1 / 

TTi + ZTTa), 7T + = -=(7Ti - MT 2 



■ ^ 

Integrating out the conjugate momenta in equation (4) and using the redefined pion 
fields, one obtains 



[dir + ] [dn_] [dn ] [da]exp 



d 4 xC 



eff 



(7) 



where £ e // is the effective Lagrangian of the system and can be written as, 



eff 



- (m 2 + ^ 2 )|tt| 2 + hdn ) 2 + Uda) 



6 



A 



-^[^ 4 + 4|7r| 4 + 7 r 4 



g^ 2 K 2 - \{m 2 + 



2 ±2 

-m <p 



24 



|2„2l 



(8) 



with |(9/7r| 2 = (dp + i//5 M o) 7r +(<9 M — 2/u5 M °)7r_ and |7r| 2 = 7r + 7r_. It should be noted that 
the chemical potential \i has entered the lagrangian in |<9/7r| 2 . The following calculations 
and discussions will be based on this effective Lagrangian. 



3. Gap equations and thermodynamic potential 

Within the imaginary time formalism of finite temperature field theory, we have 

where f3 is the inverse temperature, (3 = 1/T; the integration over the time component 
ko has been replaced by a summation over discrete frequencies. For boson there are 
u n = 2imT and n = 0, ±1, ±2, • • •. For the sake of simplicity in what follows, a 
shorthand notation is used to denote the integration and the summation. 

According to the Lagrangian (8), the tree level inverse propagators of a, 7Tq and tt± 
can be written respectively as 

D; 1 = ^ 2 + k 2 + m 2 + ^ 2 , (10) 
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D~ 



< + k 2 + m z + r 

O 



A 



{uo n + z/i) 2 + k 2 + m 2 + -0 2 . 

o 



5 

(11) 
(12) 



The equation (12) represents the inverse propagator of tt + and 7r_. As the summation 
in equation (9) is symmetric over n from — oo to +oo, uj n + ifj, and u n — ip, are equivalent 
in describing the propagators of tt + and 7r_. Here we just choose uj n + i\i to describe 
both 7r + and 7r_, and write the two propagators in the same form as in equation (12). 
According to jTHI, we can write down the corresponding effective potential at finite 
temperature as 



V(</> ,G) = im 2 2 + A 4 + 1 I lnG -i + D -i Ga 



P 



1 

+ 2 



D^Gq - 1 



2 

lntr 1 







2 

D^G 



2J0 



r lnG n 



1 



(13) 



where Go-, Go and G are the full propagators of a, 7i and ir± respectively. They are 
determined by the stationary condition (1). The chemical potential is included in the 
full propagator. ^(0, G) represents the infinite sum of the two particle irreducible 
vacuum graphs. However, in CJT formalism at the Hartree approximation we need 
only to calculate the "double bubble" diagrams with four-point vertex and treat each 
loop line as the full propagator [TUl EH] . Therefore, V 2 can be written as 

A r r I 2 A r r I 2 A r r i 2 

V 2 ( 





2 A 




2 A 




[ G a 


+ 3 


\S G 


+ 8 


/ G 


h - 




Up _ 




U0 J 



+ G a ( G + 4/ G a \ Go + ^/ Gf G . 

o J0 J 3 12 J a J0 o J0 J0 



(14) 



By minimizing the effective potential V(4>, G) with respect to the full propagators, 
we get the following set of nonlinear gap equations, 



g: 



G 1 



d: 



Do 1 



G a 



A 
2 



G 







[Go + ±[G a + ±[< 
J0 o J0 6 J0 

A 



G-i = D-i+ 2 4 /g+^/g ct + ^/Go. 

O J/3 J/3 O J/3 

Furthermore we take the similar ansatz for the full propagators as in [IT)] . 



A 

6 7/3 
A 

3 
A 



Go, 
G. 



(15) 
(16) 
(17) 



g: 



k 2 + M 2 , 



G 1 



(18) 
(19) 
(20) 



u 2 n + k 2 + Mq, 
G" 1 = K + ipf + k 2 + M 2 , 

where the effective masses M a , M and M have been introduced for a, tt and 7T± 
respectively. Substituting them into the above gap equations, we obtain a set of effective 
mass gap equations, 



M 2 



A 



A 



m + + 7T / G a + - 



A 



2 7/3 



3 7/3 



u J/3 



Gr 



(21) 
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A 
A 



A 
2 

2A 



m z + -d> z + — G+- 



X 

6 Jg 
A 



JB 



6 3 Jb 6 J /3 
Accordingly the effective potential can be written as 

1 



6 J a 



G. 

G . 



6 

(22) 
(23) 



V(<j>,M) 



-m 2 6 2 + — < 
2 y 24 



+ 



2Jp 



lnG: 



(Mi 

- [ (M 2 - m 
JB 



6 



l )G c 



InG' 



B 



h^G + V^M). 
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By minimizing the potential with respect to the order parameter 
equation, 



m 2 6 + -<p 6 + — 
r 6^ 2 



G (j 



A0 

IT 



Go 



A0 

T 



G = 0. 







(24) 

we obtain one more 
(25) 



From equations (21) — (23) and (25), the effective masses and order parameter can 
be solved self-consistently at given temperature and chemical potential. There are 
some discussions concerning the renormalization on the CJT formalism in different 
models ^UjnBlEI]- Recent investigation about the renormalization of the O(N) linear 
sigma model has been addressed in jTTj. In our discussion, the divergent parts of the 
integrals are neglected as is done in [TU]. And we take the finite parts of the integrals 
as follows: 



/ InG" 1 = / ]n[(u n + %iif + k 2 
Jb jb 



M 



1 



d 3 k 

(2tt) 3 

G = 
d 3 k 



lnfl 



ln(l 



3 - / 8(o;-/«)> 



(26) 



(2tt) 3 2a; 



B (iO n + 2/i) 2 + k 2 + M 2 

1 1 

+ 



(27) 



where u = Vk 2 + M 2 . The integration over Go and G a will be dealt with in the same 
way. With the chiral limit (e = 0), the corresponding coupling constant and the negative 
mass parameter can be taken as A ~ 125 and — m 2 w 1.8 x 10 5 MeV 2 . 

In the discussion of the thermodynamic system, the effective potential V is 
equivalent to the thermodynamic potential Q. Thus we get 



n 



1 2/2 , ^ jA 

-m q> H do 

2 ^ 24^ 



+ T 

A n. 



_fk_ 

(2^)3 
d 3 k 1 



ln(l 



-Bu>o 
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d 3 k 

(2^)3 
d 3 k 
(2^)3 



ln(l 



2 ^ j.2> 
m - -0 





d 3 k 1 



(27r) 3 cj e^o 



ln(l 



ln(l 



- (M 2 - m 2 - 



6 



d 3 k 1 

(2vr) 3 2^ 



1 



+ 



1 



/3(w+At) _ X e /3(^-M) — 1 



+ ^2, (28) 



where f2 2 = V 2 , clo- 
the relation 

P = - 



k 2 + M> 



<9p 



k 2 + M 2 and = Vk 2 + M 2 . According to 

(29) 



and equation (23), we can get the net charge density 



P 



d 3 k 



1 



This expression of density seems very similar to that of ideal gas, but actually they are 
different. Because here the effective mass M in u is a function of temperature and will 
be determined self-consistently by the gap equations. 



4. BEC and chiral phase transition 

Now we are in a position to discuss BEC and the chiral phase transition. We will study 
BEC first. It is known that when \i = M, BEC occurs. The equation (30) should be 
written as 



p = Po + P*(/3,/i = M), 



p*(/3,/i = M) 



d 3 k 



1 



(31) 



(2vr) 3 [ e P{"~M) _ 1 e /3(u+M) _ I 

where po represents the charge density of the zero-momentum state [TTJ 120]. Let us 
determine the point at which p = M is reached. From equation (30), when p is fixed, 
by solving the gap equations (21) — (23) and equation (25), we find both p and M are 
functions of T. It can be plotted in figure 1. We can see that with T decreasing, M 
decreases first and then increases, while p keeps increasing quickly and approaches M. 
At certain temperature, p catches up with M and BEC happens. From equation (31), 
we know the critical temperature T c is determined implicitly by the equation 

p = p*{(3 c ^ = M). (32) 

When T < T c , the system goes into BEC phase. The equation (31) will be solved 
together with the gap equations with p fixed. We find p and M are still functions of T 
and p(T) = M(T). They both decrease with temperature decreasing, which is indicated 
in figure 1. 

From equation (31) which is solved together with the gap equations at fixed p, we 
can get po as a function of T. This is shown in figure 2. It is clear that when p is fixed, 
at T > T c , po = 0; at T < T c , p starts to increase quickly with temperature decreasing, 
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Figure 1. p and M as functions of T at the fixed total charge density p — 0.06/m 
BEC happens at T c = 80MeV. 
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Figure 2. The ratio po/p as a function of the ratio T/T c at p = 0.06/m 3 and 
T c = 80MeV. 

which means large amount of the charged pions will reside to the zero momentum state 
and form BEC. 

For different fixed total density p, there will be different values of critical 
temperature T c and chemical potential p(T c ), so one can study the \i — T phase diagram 
of BEC. In solving the self-consistent equations, we observe that T c and p(T c ) depend 
on the order parameter 0. In low temperature regime, the order parameter is nonzero 
and determined by equation (25). When it is solved together with the gap equations 
(21) — (23) and equation (32), the critical temperature of BEC at a given density can 
be determined. The \x — T phase diagram for BEC can be given accordingly as shown 
in figure 3. The dividing line which separates BEC and the normal phase at 7^ is 
curve AC. At T > T c (the temperature at the point C), there is no solution with 7^ 
for the onset of BEC. Because in the high temperature regime the chiral symmetry is 
restored and = 0. At = the phases dividing line is curve DF as shown in figure 
3. We can see that BEC happens in chiral symmetry broken state and restored state at 
low temperature and high temperature respectively. 

In the regime of Tb < T < Tc, it seems that BEC can happen at both 7^ and 
= 0. However, if we calculate thermodynamic potential from equation (28), we find 
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under the same temperature the thermodynamic potential of ^ has lower values 
than that of <ft — 0, which means at T# < T < Tc, BEC would happen along the curve 
BC. Then the final // — T phase diagram of BEC is determined and shown in figure 4. 
In low temperature regime BEC phase and the normal phase are divided by curve AC, 
while in high temperature regime, the phases are divided by curve EF. The whole phase 
plane has been divided into two phase areas. The area above the curve belongs to the 
BEC phase, while the area below the curve represents the normal phase. 




Figure 3. The phase diagram of /z versus T for BEC. (AC and DF are the dividing 
lines which separate BEC and normal phases in low temperature regime and high 
temperature regime respectively.) 




Figure 4. The y~ — T phase diagram for BEC with no ambiguity. 



Now we will discuss the chiral phase transition in this system. By using the 
linear sigma model, there are many discussions on the chiral phase transition at 
finite temperature fUl El U2 ■ As we have considered finite pion density as well as 
finite temperature, the calculation becomes more involved. However, the chiral phase 
transition is still characterized by the order parameter 0. In our discussion, is not 
only temperature dependent but also chemical potential dependent. In fH|, the chiral 
phase transition has been discussed at finite temperature with the linear sigma model. 
By using CJT resummation scheme in Hartree approximation, the transition is found 
to be first order in the chiral limit. If we take // = 0, we can reproduce their result. 




Figure 5. as functions of <j>. Solid lines stand for T = WOMeV, and the chemical 
potentials for different lines are 260MeV, 237 MeV and 220MeV respectively (from 
top to bottom). Dashed lines stand for T — 180Mey, and the chemical potentials are 
UOMeV, lOOMeV and 60MeV respectively (also from top to bottom). 



In the case of fi ^ 0, we find the chiral phase transition is still first order. This can 
be seen if we take a look at Q as a function of 0. As the gap equations (21) — (23) 
can be solved at certain T and /i, from equation (28), the thermodynamical potential 
Q can be determined. The dependence of Q can be plotted as shown in figure 5. 
In the transition, it can be seen that there are two minima of Q which correspond to 
the stable and meta-stable states at certain T and /i. If T is fixed, by increasing /i, 
the lower minimum of 7^ 0, which stays in a stable state, will be lifted up. When it 
becomes higher than the = minimum, it changes to a meta-stable state from the 
stable state. In other words, the = minimum becomes a stable state which means 
chiral symmetry is restored. It is clear that the chiral phase transition is first order. The 
transition is usually defined at the point where the two minima become degenerate. In 
figure 5, when T = lQOMeV and T = 180MeV, we can see that the transitions happen 
at jj, — 237MeV and /j, = lOOMeV respectively. According to this way, the fx — T phase 
diagram for the chiral phase transition can be plotted as shown in figure 6. 

350 , , , , , , , 

300 ! 




130 140 150 160 170 180 190 
T (MeV) 



Figure 6. The phase diagram of /i versus T for the chiral phase transition. When 
T > lAAMeV, the area above the curve is the restored phase; the area below the curve 
is the broken phase. 
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From figure 6, we can see that the transition temperature is 184MeV for \i = 0, 
which was obtained in ^U]. Here when \i ^ 0, the transition temperature will decrease 
with the chemical potential increasing. At certain temperature (T m 144MeV), the 
phase curve of the chiral phase transition will terminate, because at this time fi = M 
and BEC will happen. As a result the chiral phase transition can not take place further 
at lower temperatures. When T > 144MeV, the area above the curve is the chiral 
symmetry restored phase. However, the phase diagram is not complete, because we do 
not know what the phase is when T < 144MeV . Now let us recall the \i — T phase 
diagram of BEC (figure 4). If we combine the two phase diagrams (figure 4 and figure 
6) into one figure, we can see the different phase areas clearly. This is shown in figure 
7. When the chiral phase transition terminates, it means that the system will go into 
BEC phase. Furthermore, the normal phase area in figure 4 is divided into the chiral 
symmetry broken phase and the chiral symmetry restored phase. When the system is 
in a chiral symmetry broken phase, if we decrease temperature with density fixed, BEC 
will occur at ceratin temperature; if we increase temperature, the system may go into 
a chiral symmetry restored phase. In the restored phase, at relatively low densities, if 
we decrease temperature, the chiral phase transition may take place; at relatively high 
densities, if we decrease temperature with density fixed, BEC will take place. 




25 50 75 100 125 150 175 200 
T (MeV) 



Figure 7. The combined phase diagram of \i versus T for BEC and the chiral phase 
transition. The CSB and CSR stand for chiral symmetry broken phase and chiral 
symmetry restored phase respectively. 

We noted that in |Hj, BEC and the spontaneous symmetry breaking in interacting 
Bose system have been discussed with the O(N) model. In their discussion, a chemical 
potential /i associated with 0(2) symmetry is introduced when N is even. There are N/2 
such chemical potentials for N scalars. For only one chemical potential and the large-N 
limit, their result shows that the transition to broken symmetry phase is just the onset of 
BEC in the the ultra-relativistic regime. This is because the scalar field associated with 
\x acquires a nonzero vacuum expectation value at the point that \l — M is reached when 
temperature is lowered. Thus the symmetry, which is actually the 0(2) symmetry, is 
broken at this time. In our discussion, we have studied the chiral symmetry restoration 
and BEC in the 0(4) linear sigma model. In this model, the four scalar fields have real 
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physical meaning, which are three pion fields and one sigma field. The chiral symmetry 
breaking in the vacuum is chosen along the a direction. The chiral phase transition is 
characterized by the order parameter <fi. The chemical potential here is associated with 
the pion fields. The BEC means the BEC of the charged pions. The relation between 
BEC and the chiral phase transition can be illustrated by the \i — T phase diagram as 
indicated in figure 7. 

5. Conclusion 

In summary, we have discussed Bose-Einstein condensation and the chiral phase 
transition for the linear sigma model in the chiral limit. The effective potential of 
the linear sigma model is calculated within the CJT resummation scheme in the 
Hartree approximation at finite temperature and finite pion density. By solving the 
self-consistent equations we have discussed BEC in the interacting pion system. The 
fx — T phase diagram of BEC is presented. At relatively low temperature, BEC happens 
in chiral symmetry broken state, while at relatively high temperature, BEC happens 
in chiral symmetry restored state. The phase plane has been divided into the BEC 
phase and the normal phase. Then we have discussed the chiral phase transition and 
also plotted the \x — T phase digram of the chiral phase transition. By comparing it 
with the BEC phase diagram, we find the chiral phase transition happens in the normal 
phase and further divided the normal phase into the chiral symmetry broken phase 
and restored phase. So the whole phase plane has been divided into three phase areas: 
the BEC phase area, the chiral symmetry broken phase area and the chiral symmetry 
restored phase area. BEC can happen either from the chiral symmetry broken phase or 
from the chiral symmetry restored phase. When the system goes into BEC phase, there 
is no chance for the onset of the chiral phase transition. 
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